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Abstract 



Near-extreme Reissner-Nordstrom-anti-de Sitter black holes are unstable against the condensation 
of an uncharged scalar field with mass close to the Brcitenlohner-Freedman bound. It is shown that a 
similar instability afflicts near-extreme large rotating AdS black holes, and near-extreme hyperbolic 
Schwarzschild-AdS black holes. The resulting nonlinear hairy black hole solutions are determined 
numerically. Some stability results for (possibly charged) scalar fields in black hole backgrounds 
fH are proved. For most of the extreme black holes we consider, these demonstrate stability if the 

"effective mass" respects the near-horizon BF bound. Small spherical Reissner-Nordstrom-AdS 
black holes are an interesting exception to this result. 
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1 Introduction 

Charged black holes in anti-de Sitter space have played an important role in recent discussions of 
superconductivity based on the gauge/gravity correspondence. The simplest model of this phenomenon 
consists of gravity coupled to a Maxwell field and a charged scalar field [T] . The non-superconducting 
phase is dual to the Reissner-Nordstrom-AdS black hole solution, which is stable at sufficiently high 
temperature. However, at temperature below a certain critical value, this solution is unstable against 
condensation of the scalar field. At these low temperatures there exists a new charged black hole 
solution, with a non-trivial scalar field present. This describes the superconducting phase. 

Surprisingly, the condensation can occur even when the scalar field is uncharged [2j. A nice way 
of understanding this is to consider the special case of an extreme (planar) Reissner-Nordstrom-AdS 
black hole. This has AdS2 x R^ near-horizon geometry. The Breitenlohner-Preedman bound [H H] 
required for stability of the AdS2 lies above that required for stability of the asymptotic AdS^. This 
suggests that an uncharged scalar field which satisfies the AdS^ BF bound but violates the AdS2 BF 
bound will be unstable in the extreme RN-AdS backgroundjj This is confirmed by solving numerically 
the scalar equation of motion in the black hole background [U E] . 

The motivation for the present paper is the observation that extreme, rotating, AdS black holes 
also have near-horizon geometries containing an AdS2 factor. Hence one might expect that a scalar 
condensation instability will afflict these black holes toojj We shall confirm that this is indeed the 
case for black holes sufficiently close to extremality. 

We shall consider the (topologically spherical) rotating AdS^ black hole solutions discovered by 
Hawking, Hunter and Taylor (HHT) ^ (sometimes called Myers-Perry-AdS). In general these have 
two angular momenta but we shall set these equal since the solution then depends non-trivially only 
on a single radial coordinate r. These solutions can be parameterized by the horizon radius r_|_ and the 
angular velocity CIh- At extremality, the AdS2 BF bound is above the AdS^ bound if r+/£ > 1 (where 
£ is the ^^55 radius), i.e., for "large" black holes. In this case, we might expect scalar condensation 
to occur if the scalar mass /j, lies between the two BF bounds. 

We have determined numerically the values of r+ , 0,h and the scalar mass /i for which there exists 
a time-independent solution of the scalar field equation of motion which preserves the symmetries of 
the background. Such a solution is expected to arise at the threshold of the instability. Figure [l] 
summarizes the result of this calculation. The condensate appears only when the black hole is very 
close to extremality: typical temperatures are Th£ ~ 10^'^. 



^Ref. [2] attributed this argument to M. Roberts. 

^It was shown in Ref. [6] that the charged scalar condensation instabihty of static charged AdS black holes extends 
also to rotating charged AdS black holes. In our work, neither the scalar field nor the black hole is charged. 
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Figure 1: Scalar condensation in the HHT black hole with equal angular momenta. We consider a 
free scalar field whose CFT dual is an operator of dimension A. The vertical axis is the difference of 
Qh from its extreme value 5Q = 17^* (r^-) — Qjj- The surface corresponds to values of ^h, ^+ and 
A for which there exists a regular, time-independent solution of the scalar equation of motion. The 
surface extends to arbitrarily large r_|_. The region enclosed by this surface corresponds to values of 
Qh, ?*+ and A for which the black hole is unstable against scalar condensation. The intersection of 
the surface with the 60, = plane corresponds to saturation of the near-horizon AdS2 BF bound for 
extreme black holes. The black curve is for A = 2, i.e., saturation of the AdS^ BF bound. 



The analogy with the charged case suggests that the endpoint of the instability should be a 2- 
parameter rotating black hole solution with scalar hair. This new solution should have the same 
symmetries as the HHT solution, so it will have equal angular momenta. We have determined certain 
1-parameter subfamilies of this solution numerically (taking fi to saturate the AdS^ BF bound). We 
find that the hairy black holes have higher entropy than the HHT black holes with the same mass 
and angular momenta. Hence it is consistent with the second law for the HHT black hole to evolve to 
our new hairy solutions. Our numerical results suggest that the entropy of the hairy black holes will 
vanish as the temperature approaches zero, just as in the charged case [H iQtITO]. 

We have extended the linear analysis to the four-dimensional Kerr-AdS solution. Again, we found 
that a scalar condensation instability occurs in a region close to extremality. Since this solution 
depends non-trivially on two coordinates (r and 0), we did not attempt to construct numerically the 
corresponding hairy black holes. Note that it has been proved that a scalar field that respects the 
AdS^ BF bound is stable in a sufficiently slowly-rotating Kerr-AdS background [H]. Our results show 
that this conclusion does not extend to rapidly rotating Kerr-AdS black holes. 

Our results have implications for (conjectured) no hair theorems. Consider AdS gravity coupled 
to an uncharged scalar field. Ref. [12] conjectured that static spherically symmetric hairy black hole 
solutions do not exist if the theory admits a positive energy theorem with boundary conditions that 
preserve the symmetries of AdS. Our results demonstrate that this conjecture cannot be extended to 
non-static black holes. 

For Reissner-Nordstrom-AdS, the instability against condensation of an uncharged scalar can co- 
exist with an instability against condensation of a charged scalar. In the rotating case, the analogue 
of the scalar field charge is the angular momentum quantum number m of the scalar field, which we 
have set to zero. However, for m ^ 0, black holes with ^h^ > 1 exhibit a superradiant instability 
|13j . This is the analogue of charged scalar condensation. Now, since extreme black holes always have 



^H'^ > 1, and scalar condensation with m = occurs only very close to extremality, it follows that the 
scalar condensation instability studied in this paper always coexists with the superradiant instability. 
Moreover, the new hairy black holes that we find have Qh^ > 1- Therefore they are likely to suffer 
from a superradiant instability themselves. The endpoint of this kind of instability is unknown. 

Two other topics are studied in this paper. First, we note that a nice toy model for the study of 
uncharged scalar field condensation is the Schwarzschild-AdS black hole with hyperbolic spatial slices. 
This solution has a regular extreme limit (in contrast to the spherical or planar Schwarzschild-AdS 
black hole) with an AdS2 in the near-horizon geometry. We determine the range of values for r+ and 
fj, for which scalar condensation occurs. The resulting 1-parameter nonlinear hairy black hole solution 
is determined numerically (for particular /u)r| 

Second, we shall present some stability results concerning scalar condensation. The idea is to 
construct an "energy" functional for the scalar field (which might not coincide with the usual energy) 
which is manifestly non-increasing in time. If one can show that this energy is positive definite 
then an initially small fluctuation in the field must remain small, i.e. condensation does not occur. 
Demonstrating positivity is non-trivial because we are mainly interested in the situation in which 
IJ? < 0. Fortunately, this problem has been addressed in Ref. |llj for spherical Schwarzschild-AdS 
black holes so we can generalize this method. For the extreme black holes considered in this paper, we 
show that the energy is positive definite if, and only if, the AdS2 BF bound is satisfied. Our numerical 
results demonstrate that an instability appears as soon as the bound is violated. For non-extreme 
black holes, our method gives a lower bound on /i^ that guarantees stability although this bound is 
not sharp. 

The method can also be applied to a charged scalar field in a Reissner-Nordstrom-AdS background. 
For planar, hyperbolic or sufficiently large spherical black holes, our energy functional is positive 
definite at extremality if, and only if, the near-horizon "effective mass" satisfies the AdS2 BF bound. 
For the planar case, Ref. [5] found numerically that an instability occurs precisely when this bound 
is violated. Small spherical extreme RN-AdS black holes are an interesting exception: for these our 
stability bound is more restrictive than the near-horizon BF bound and indeed Ref. [9] found an 
instability even when the near-horizon BF bound is respected. 

This paper is organized as follows. We state the problem in Section [2j We then discuss in Section [3] 
the hyperbolic Schwarzschild-AdS black hole as a toy model to illustrate the scalar condensation 
phenomenon. Next we discuss rotating AdS black holes: cohomogeneity-1 solutions in Section [4] and 
also the four-dimensional Kerr- AdS in Section [5} In Section |6] we present our stability results for scalar 
condensation. 

2 Scalar condensation 

Consider a scalar field with mass n in an asymptotically AdS^ background. It satisfies the Klein- 
Gordon equation, 

V2$-^2^ = 0, (1) 

and always decays at infinity as [31 ^ 



4(+) /4(-) d-1 lid- 1)2 

Hr)^-^ + -^, where A± = ^- ± J^— ^ + ^^^^ (2) 
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and I is the AdS radius. Stability of the AdS background requires reality of A-t, i.e. that the mass of 
the scalar field must obey the Breitenlohner-Freedman (BF) bound [3l H] 



(d-1) 



^ >/i \bp = -Tj^. (3) 



4^2 



^We note that previous work has found a similar solution analytically by including a suitably chosen scalar potential 
|14| . We emphasize that we are working with a free (and minimally coupled) scalar. Another difference is that we use 
AdS invariant boundary conditions at infinity, which does not appear to be the case for the solution of Ref. [14] . 



For scalars with mass above this BF bound and below the unitarity bound, i.e. 

^''\BF<^''<^'\n^t^ with ^\^^^^JAZ^ + }.^ (4) 

there is a choice of boundary conditions: one can impose either that <& decays as r~ + or as r~ , 
since both give normahzable solutions. In the AdS/CFT correspondence, this choice dictates whether 
the operator dual to ^ has dimension A_(_ or A_ [15j. For masses above the unitarity bound only the 
mode with the faster fall-off, r~ +, is normahzable. 

Now consider an extreme, asymptotically AdSd black hole whose near-horizon geometry contains 
an AdS2 factor with radius R. The BF bound associated to this AdS2 is 



Hence if 



2^21 _ '- /K\ 

^ - ^ iNHBF - 4^2- W 
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^ \nhbf>^' ^/^ \bf (6) 

then the asymptotic AdSd will be stable but the near-horizon geometry is unstable. This suggests that 
the full black hole solution will be unstable against scalar condensation |2|. This has been confirmed 
by numerical calculation in certain cases |5]. 

3 Toy model: hyperbolic Schwarzschild-AdS 

3.1 Introduction 

The simplest black hole that exhibits the scalar condensation phenomenon is the Schwarzschild-AdS 
solution with hyperbolic spatial sections. Consider the general Schwarzschild-AdS solution in d di- 
mensions: 

ds^ = -f{r)de + f{r)-^dr^ + r^dT^, (7) 

where (iS| is the metric on a unit sphere {k = 1), hyperboloid {k = —1), or flat space (A: = 0) that 
satisfies Rab = k{d — 3)gab- The function /(r) is given by 

(^d-3 _ ^rf-3) ^d-i _ ^rf-i 

where the horizon is at r = r+. 

For k = 0,1, the horizon is always non-degenerate. However, for k = —1 a black hole solution 
exists for 

r+>r''f, with rf = y?H|^, (9) 

and has Hawking temperature 

{d - l)rl - {d - 3)f , , 

An extreme (zero temperature) regular solution is present when the bound in (lo]) is saturated. This 
extreme solution has near- horizon geometry AdS2 x H'^~'^, where the AdS2 has radius R = ll\/ d — 1. 
Hence scalar condensation seems likely when the scalar mass lies between the AdSd and AdS2 BF 
bounds: 

If n lies strictly above the AdSd BF bound then the stability argument that we shall present in 
Section [6] below demonstrates that the scalar field is stable for sufficiently large r_|_. However, if /x also 
lies below the near-horizon BF bound for the extreme black hole then we expect the scalar field to 
become unstable as extremality is approached, i.e. as r+ decreases. At the value of r+ corresponding 



to the threshold of instabihty, we expect there to exist a time-independent solution <l>(r) of the scalar 
equation of motior|j 

<^"ir) + (^ + 7^) <^'ir) - f^fir) = • (12) 

We shall confirm the existence of this solution numerically. 

We note that the entropy and energy of the k = —1 black hole are given by 

where Vs is the volume of the hyperboloid (which we shall assume to be compactified) . 

The solution with r_|_/^ = 1 is special: this solution is locally isometric to AdS^ and has zero 
energy. Solutions with r+/^ < 1 have negative energy with the extreme solution having the lowest 
energy. 

3.2 Numerical results: linear 

For definiteness, we set d = 5. For scalar mass below the unitarity bound we shall consider both 
choices of boundary condition. For this reason, it is more convenient to work with the dimension A 
of the operator dual to the scalar field than with /x. For given A we wish to determine the value of 
r+ for which there exists a solution $ = $(r) that is regular at the horizon and decays as r~ at 
infinity. We do this using a standard shooting method where we construct power series solutions near 
the horizon and at infinity up to seventh order and then evolve the solution away from these points 
numerically. Fixing the amplitude of the solution at spatial infinity, we impose matching conditions, 
namely the continuity of ^ and <!>', at some fixed value of r > r+. We find that matching is possible 
only for a particular value of r+/£. 

The results are presented in Figure 2 The black hole is extreme for r-|_ = r'j!^* = 2^^'^ ~ 0.707. 
The near-horizon BF bound (which lies above the unitarity bound) corresponds to A = 2+\/3 ~ 3.732. 
This is indeed the limiting value of the curve at the extremal blue point jj Hence our results agree 
with the expectation that scalar condensation occurs for the extreme solution if, and only if, the 
near-horizon BF bound is violated. For A > 2 + v3, we find no solution, so there is no scalar 
condensation instability for any r_|_. However, as A decreases through 2 + v3, first the extreme 
black hole becomes unstable and then the critical temperature below which the black hole is unstable 
increases monotonically. 

Saturation of the AdS^ BF bound corresponds to A = 2. Our numerical results suggest that the 
threshold of instability for A = 2 occurs at r_|_ = L One can confirm analytically that the black 
hole with r_|_ = £ admits the regular solution $(r) = £^/r^ if A = 2. As mentioned above, the 
r_|_ = £ solution is special: it is locally isometric to AdS^. Our results imply that this solution is 
unstable against scalar perturbations with A < 2. But AdS^ is stable! The point is that our solution 
is isometric only to part of AdS^. If one analytically extended an unstable mode of our solution to 
global AdS^ then either it would not be regular everywhere or it would not satisfy the appropriate 
boundary conditions at infinity. Hence it would not correspond to an instability of AdS^ so there is 
no contradiction. 

3.3 Numerical results: nonlinear 

In this subsection, we construct the family of hairy black holes that branches-off from the original 
Schwarzchild-AdS black hole family in a phase diagram of static and neutral solutions with an hyper- 
bolic horizon topology. 



^ See Refs. [161 117| for previous studies of scalar fields in the hyperbolic Schwarzschild-AdS spacetime. 

^ This blue point is a limit of the curve, it does not belong to the curve. This is because, at extremality, there is no 
regular time-independent solution associated to the threshold of instability. The scalar field diverges at the horizon if 
one takes a limit of the time-independent solution associated to the threshold of instability of non-extreme solutions. 




Figure 2: Threshold of the scalar condensation instability in the d = 5 hyperbolic Schwarzchild-AdS 
black hole. The plot shows the critical dimensionless horizon radius r^/£ where the instability sets in 
as a function of the conformal dimension A of the CFT operator dual to the scalar field. Points below 
the curve correspond to instability. The blue dot corresponds to the extreme black hole. The critical 
value of A at this point is precisely the value predicted from the near-horizon BF bound. The black 
dot corresponds to the AdS^ BF bound, i.e., A = 2. 



We employ two different methods that agree in the regime where they have to. First, we follow 
an analytical perturbative approach. Starting with an appropriate ansatz for the hairy black hole 
solution, we expand out the gravitational and scalar fields in a power series of a small parameter e. 
This is a measure of the value of the scalar field at infinity and indicates the deviation of the hairy 
solution from the unstable Schwarzchild-AdS black hole. We do this perturbative expansion up to 
15**^ order. To obtain the exact result, we then do a full non-linear numerical construction of the hairy 



black hole solution. We present our results in subsection 3.3.3 where we describe our findings in a 
phase diagram of static and neutral solutions with an hyperbolic horizon topology. 



3.3.1 Hairy black hole ansatz. Equations of motion 

We want to look for hairy black hole solutions with an hyperbolic horizon topology. As in the linear 
analysis, we consider d = 5. Therefore we take the following ansatz for the gravitational and scalar 
fields, 



ds"^ = -/(r)e^(^) dt^ + f{r)-^dr^ + r^d^l^ 
^ = (I)(r) . 



(14) 



When $(r) = 0, xi^") = and /(r) is given by ([8| with k = —1, the system describes the hyperbolic 
Schwarzchild-AdS black hole. This ansatz solves the equations of motion derived from the Einstein- 
scalar action, 
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if the following equations are satisfied 
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^2^2 



$"(r) + 



[(6(2r2 - f) + 3ff{r) - r^n^f^{rf) $'(r) - Srfi'^f^r)] = , 

2 4r 



3Prf{r) 
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(15) 

(16a) 
(16b) 



r . 



X'{r) = -^'{ry. (16c) 



Given <^(r), x(^) can be easily found through the integration of (16c). So the problem at hand reduces 
to determining <&(r) and /(r). 

The boundary conditions to be imposed are as follows. We use the scaling symmetry of the 
equations of motion, e^ — )• X'^e^ and t — )• At (with the other variables/fields unchanged) to set x = 
at the asymptotic boundary without loss of generality. At the horizon we demand that x is regular. 
The function /(r) must vanish at the horizon, and we take this condition as our definition for the 
location of the black hole horizon. We require the hairy black hole to be asymptotically AdS and 
thus /(r) must behave as r^/£^ — 1 at large distances. On the other hand, the scalar field has to 
be regular at the horizon. Its asymptotic boundary condition is determined by the requirement of 
normalisability at infinity. In our linear analysis we found that a hairy black hole should exist only for 
scalar masses above the BF bound and below the NH BF bound. Therefore we restrict our attention 
to this range of masses. For concreteness we consider the simplest case fi = hIbf (and thus A = 2) for 
which an analytical perturbative construction is possible and where the numerics simplify considerably. 
Summarizing, the boundary conditions for the non-linear problem are: 



r 



2 



-4\ 



/ =0, /K ^1^-1; X =0(1), xK ^o{r 

$1 =0(1), $1 _^ ^(/-or-^ (17) 

lr=r-|- ^ ' ' Ir— >cxD ^" ^ ' 

For an explicit construction of the hairy black hole we should still select a particular hyperbolic 
Schwarzchild-AdS black hole to which the hairy black hole should reduce in the limit where the scalar 
field vanishes, i.e. a particular value of r_|_ in ([s]) (with k = —1). We are taking A = 2 so will find 
the hairy black hole that bifurcates from the hyperbolic Schwarzchild-AdS black hole with r^ = i. 



This choice is not arbitrary. Indeed in this case we can solve analytically the equation (12) that gives 
the perturbative leading order scalar field. This can then be used to generate the next-to-leading 
order expansion terms. Ultimately, we will then use this perturbative result as the seed solution in 
the relaxation method employed to construct numerically the exact hairy black hole. To sum up, our 
selection of parameters is 

3.3.2 Perturbative and numerical construction of solution 

To construct perturbatively the hairy black hole, we expand all the unknown functions of the system, 
namely r_|_, /(r) and <I>(r) in a power series of e, 

r+(e) = iJ2P^^ '"' ' /(^' ^) = E /2i(0 e'^ , Hr, e) =^2 ^2,+i(r) e'^+' , (19) 

i=o j=o j=o 

where e is a measure of the scalar field at infinity. Using a standard perturbation theory strategy. 



we plug these expansions in the equations of motion (16) to solve for the coefficients p2j, f2j{f) and 
$2j+i(^)- In short, the procedure at each order is the following. Equation (16b) is used to solve for 
/2j(r), up to an integration constant. The boundary condition /(r^) = must be satisfied at each 
order and is used to find p2j , which defines the horizon location where some of the boundary conditions 



will be imposed. We can now use (16a) to solve for (^2j+i{f), again up to integration constants. The 



boundary conditions (17) are now imposed to constraint the integration constants generated in the 



process. Finally, x(^) also has an expansion in e that can at this point be easily found by direct 



integration of (16c). 



At leading order, n = 0, the scalar field is given by ^(r) = el- /r + 0{e ) while /(r) = r /£ — 1 + 
0{e^) and r+ = l + 0{e^). That is, we are in the linear regime discussed in the previous section, where 



the back-reaction in the gravitational field is neglected and $(r) solves (12) (note that the analytical 
solution is possible only because we perturb the black hole with r_(_ = i). We will find that for our 
purposes it will be enough to do the expansion up to order n = 14. The explicit expression for the 
expansion functions {p2j ,f2j{r) ,^2j+i{r)} up to 15**^ order would be far too cumbersome and not 
illuminating at all. For illustration purposes we just give the result up to the 3^'^ order, 
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e^' + Oie^) 



(20) 



The result of the perturbative analysis up to 15*^^ order will be presented in the phase diagrams of 
subsection 13.3.31 

The most efficient way to get the full information on hairy black holes in the phase diagram 
of solutions is to resort to a full non-linear numerical approach that solves non-perturbatively the 
equations of motion (16). Our numerical strategy to find the exact hairy black hole is to use a 
standard relaxation method [18j, with a spectral discretization of the integration grid [19J. Again, we 



have to solve the system of coupled ODEs (16a) and (16b) subject to the boundary conditions (17), 



and for the selection of parameters (18). 



The implementation of the boundary conditions using the spectral discretization is simpler if we 
implement the following redefinition of the fields. 
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in which case the boundary conditions (17) reduce to 
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(22) 



=r+ 



The factor ^'^/r^ in the relation for qf in (21 ) guarantees that the field redefinition does not introduce a 
new critical (singular) point in the differential equations to be solved (i.e. besides the boundary critical 
points). Our numerical results for the hairy black hole will be presented in the next subsection. 



3.3.3 Phase diagram 

In this subsection we present the properties of the hairy black holes that we constructed using the 
perturbative and numerical approaches described in the previous subsection. We present the scalar 
condensate, the temperature and the entropy of the hairy black hole as a function of its energy. 

To compute the energy of the hairy black hole in AdS we use the Astekhar-Das formalism |20j . 
The temperature, entropy and energy of the hairy black hole are given by 
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(23) 



where the energy is measured with respect to the AdS background. Vs is the volume of the (compact- 
ified) hyperboloid of unit radius, (J)q£'^ is the 0{r^'^) coefficient in the large r expansion of <^(r) and 



Tq is the 0{r~'^) coefficient in the large r expansion of f{r). We will compare these thermodynamic 
quantities with those for the family of hyperbolic Schwarzchild-AdS black holes. Henceforth we shall 
set Vs = 1, equivalently E and S can be regarded as energy and entropy per unit hyperboloid volume. 



In Fig. ^ we plot the value of the condensate at infinity, (/>o/^^ as defined in (17), as a function 
of the dimensionless energy E/(^ of the hairy black hole (equivalently: we are plotting the vev of 
the operator dual to our scalar field). The condensate vanishes in the limit where the hairy black 
hole reduces to the Schwarzchild-AdS black hole with r_|_ = i, which has zero energy. The condensate 
then increases as the energy grows negative. The red curve describes the results obtained using the 



analytical perturbative construction of section 3.3.2 while the blue curve describes the exact numerical 



results using the relaxation method of section 3.3.2 The matching of these two curves is very good 



for absolute values of the energy that are not too large, as it should be. 
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Figure 3: Asymptotic value of the condensate, (/>o/^^ as defined in (17) (equivalently: vev of A = 2 
operator dual to the scalar field), as a function of the dimensionless energy E/(^ of the hyperbolic 
hairy black hole. The red curve is the perturbative result, the blue curve is determined numerically. 

In Fig. H we plot the dimensionless temperature, Tjjl (left), and dimensionless entropy, S/l^ 
(right), of the hairy black hole as a function of its dimensionless energy E /i"^ . Again, the red (blue) 



curve describes the results obtained using the perturbative (numerical) construction of section 3.3.2 



Once more these two curves coincide for absolute values of the energy that are small. The perturbative 
construction up to the 15**^ order in the expansion starts deviating from the numerical curve only for 
energies smaller than E/i'^ ~ —0.02 (say). In these plots, we also display in a dashed black line 
the temperature and entropy of the hyperbolic Schwarzchild-AdS black hole family. We display this 
family only for r^* < r+ < £ which corresponds to a range of (non-positive) energies where it co-exists 
with the hairy black hole. At r_|_ = £, we have the Schwarzchild-AdS black hole with zero energy, from 
which the hairy black hole solution bifurcates. As we decrease the energy of the Schwarzchild-AdS 
family, the temperature and entropy starts decreasing until the extreme black hole with r+ = r*p*, 
T}j = and non-vanishing entropy is reached. 

Our results above show that hyperbolic Schwarzschild-AdS black holes with rj^/l < 1 {E < 0) are 
unstable against condensation of a free scalar field dual to a A = 2 CFT operator. The threshold 
mode of this instability is identified as a black point in Fig. [2} In the phase diagrams of Fig. [4j this 
zero-mode signals a bifurcation point at E' = in the Schwarzchild-AdS family of solutions into a 
new branch of hairy black holes. For fixed energy where the two families co-exist, the hairy black 
holes always have higher temperature and entropy than the Schwarzschild-AdS black holes. So we 
have a second order phase transition. As the energy decreases, the temperature and entropy of the 
hairy black hole decrease. This family of solutions ends in a zero temperature state where both the 
temperature and the entropy of the hairy solution vanish. This configuration is represented by a green 
point in the phase diagrams. Using the information from both phase diagrams we find this state to 
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have an energy of E/l"^ ~ —0.02487 with an estimated error based on both plots of 10~^ %. 

It is worth commenting on the nature of the zero temperature hmit of these hairy black holes. 
The fact that the entropy appears to vanish in this limit is good evidence that the limiting geometry 
is not an extreme black hole. Further evidence for this comes from the fact that if a hairy extreme 
black hole did exist then the scalar field would be constant in its near-horizon geometry. But then the 
scalar equation of motion would imply that the scalar field must be zero in the near-horizon geometry 
(this argument is due to Ref. [21j). This would imply that the scalar must vanish at the horizon of 
the full extreme black hole. However, our numerical results indicate that the value of the scalar field 
at the horizon increases as the temperature decreases. Hence the zero temperature limit cannot be an 
extreme black hole. A similar conclusion has been obtained in the case of charged scalar condensation 
in Reissner-Nordstrom-AdS [8j. 
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Figure 4: Dimensionless temperature, Th^ (left), and dimensionless entropy, S/i'^ (right), of the 
hyperbolic hairy black hole as a function of its dimensionless energy E/i"^ . The scalar is at the BF 
bound: //^^^ = —4 (A = 2). The red curve is the perturbative result, and the blue curve is the 
numerical result (it ends in a zero temperature and entropy configuration indicated by a green dot). 
The dashed black curve corresponds to the hyperbolic Schwarzchild-AdS black hole. This terminates 
at an extreme solution with non-zero entropy. Note that the hairy black hole always has greater 
entropy than the Scharzschild-AdS solution with the same mass. 



4 Rotating black holes with equal angular momenta 

4.1 Introduction 

The HHT black hole j7] with equal angular momenta is the d = 5 case of a family of "Myers-Perry- 
AdS" black holes j22l|23j with cohomogeneity-1. These solutions exist in odd dimensions d = 2A^-|-3, 
and the metric can be written as|£] 



ds' 



fir) 
h{r) 



dt^ + 



dr'^ 



+ r'^h{r)[dtp + Aadx" - n{r)dtf + r'^ gatdx" dx'' 



(24) 



^The radial coordinate used here can be related to the standard Boyer-Lindquist radial coordinate of [221123) through 
->■ (r^ + a^) [1 — J2-] . The functions (/, h) used here are related with the functions {F, G, H) used in [13] through 



f{r)/h{r) = F{ry, f{ry = G{rf and r^/iM = H{ry 
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where 



2 ^2N / „2\ ^2N„2 ^2N „2 ^2N , 



/(0 = i + i-3i 1-^ +:5^, Mr) = i + :3&, m 



£2 ^2N \ £2 / ' ^2(Af+l) ' ^ ' ' j.2{N+l)' "^ ' ^2{N +1) j^(y.y 

and cjah is the Fubini-Study metric on CP^ with Ricci tensor Rah = 2{N + l)gab , and A = Aadx"' 
is related to the Kahler form J by dA = 2J. Surfaces of constant t and r have the geometry of a 
homogeneously squashed 5"^ ^^, written as an S^ fibred over CP . The fibre is parameterized by the 
coordinate tp, which has period 27r. 

The spacetime metric satisfies Rf^y = —£~'^{d — l)g^u- Asymptotically, the solution approaches 
anti-de Sitter (AdS) space with radius of curvature i. The event horizon is located at r = r+ (the 
largest real root of f{r)) and it is a Killing horizon of ^ = 9t + ^nd^ , where the angular velocity of 
the horizon is given by: 



^2Ar 1 / ATf2 

The solution saturating the bound in the angular velocity corresponds to an extreme black hole with 
a regular, but degenerate, horizon. Note that this upper bound in ^h^ is always greater than 1 and 
tends to the unit value in the limit of large r^/i. 

We shall find it convenient to parameterize the solution in terms of (r+,il//) instead of {rM,^) 
through the relations, 

2(Af+l) / 2 ,f2\ r2/20„ 

^A^ - ^2_^2 _ „2 (^^ + £2) ' "-^^+£2- ^^^i 

The temperature of the black hole. 



^ / ri+£2 Ar^2_(^^i),^(^2^^2_i) 

^^ y ^2 _ ^^ (f^2^^2 _ 1) 27rr+^2 ' y^n 

vanishes for the extreme configuration. The mass E and angular momentum J, defined with respect 
to 5^, are [21] (taking G = 1) 



Svr ^'' \ 2 2(r2_+£2)2y ' g^^ r^ + £2 

where A2n+\ is the area of a unit (2A^ + l)-sphere. 

Once again, we shall look for a time-independent solution of the scalar wave equation that we 
expect to signal the onset of instability. We shall assume the scalar field to be axisymmetric (with 
respect to 5^), i.e., it preserves the symmetry associated to the rotation of the background. However, 
we shall allow for the possibility that the scalar breaks the symmetry of CP . We take the field to 
have the separable form <!>« = '&(?')Yk(x). Here, Y^ is a scalar harmonic of CP^ , with Laplacian 
eigenvalue given by — 4k(k + A^), k = 0, 1, 2, . . .. The scalar wave equation reduces to 

^ (0 + [-^— + J^) ^'(0 - ^ ^2;(^) -Hr) = . (29) 



In the extreme case, the near-horizon geometry is homogeneous and takes the form of a fibration 
over AdS2- Prom the above equation we can read off an "effective mass" for the scalar field near the 

,2 



horizon: ^g^ j- = iJp + 4«;(k -|- N)/r\. The near-horizon BF bound (determined by the radius of the 



AdS2) is 



^Iff > ^'\ZL = -^ (^^ + 2 + AT ^ ) (30) 



+ 
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We are interested in whether it is possible to violate this bound whilst respecting the AdSd BF bound 
([3|. A short calculation reveals that this is possible if, and only if, 



> 



8k{k + N) 



'1 + 



N{N + l) 



(31) 



So we expect an instability only for black holes larger than the AdS length i. Furthermore, modes 
with K = should be the most unstable since they have the lowest ^J'1ff■ Written out explicitly, we 
expect an instability of the extreme black hole if the scalar mass fj, lies in the range 



„2| ^ „2 . , 2|(MP) _ Mf^ + N) 

/^ \bF ^ M ^ /^ \nHBF ^2 



(32) 



Another interesting 5d black hole solution is the supersymmetric AdS^ black hole of Ref. [25j . In this 
case, we find that the near-horizon BF bound coincides with the AdS^ bound so there is no reason to 
expect a scalar condensation instability. We shall not consider this solution further. 



4.2 Numerical results: linear 



4.2.1 Numerical methods 



We will solve Eq. ( 29 ) numerically in search of the onset of the scalar condensation. This is a boundary 



value problem. The appropriate boundary conditions can be found through a standard Frobenius 
analysis. At the horizon, regularity of the solution requires that the scalar field must be finite, while 
at the asymptotic region the scalar field decays as 



$(r) ~ $*^+) 



1 



^2A+(A+-2iV) 
4r2 (A+ - N) 



+ 



$(-) 



1 



'2A_(A_ -2N) 
4r2_(A„ - N) 



(33) 



Stability of the AdS background requires that the scalar field must obey the Breitenlohner-Freedman 
bound ([3|. Both fall-offs in (33) can be normalizable so we can impose one of the two boundary 
conditions: <!>(+) = or <!>(") = 0. 



At this point we make the important observation that the radial equation (29) can be written as 
a generalized eigenvalue problem in i7^. That is, we can cast (29) in the form 



L{r)^{r) = f^J^r A(r)$(r), 



(34) 



where L{r) and A(r) are both second order differential operators that do not depend on VLh- For fixed 
r+ and A, our strategy for finding a zero-mode for the scalar condensation will be to determine the 
eigenvalue ri//(r+)2 for which there exists a solution of (34) that is regular at r = r+ and satisfies the 
appropriate boundary condition (either $(+) = or <!>("■' = 0) at r = oo. This strategy is motivated 

More 



by the availability of numerical techniques for solving eigenvalue equations of the form (34) 
concretely, we solve this equation using a spectral numerical method. This method already proved 
to be suitable to study perturbations in black hole backgrounds [26^, "ST", ^28] . The application of the 
method is simpler for Dirichlet and/or Neumann boundary conditions. Therefore we introduce the 
following perturbation functions. 



,(±)/ 



$(r) 1 



^ 4r^(A±-iV) ' ' 



\ r / 



r 



(35) 



and we impose the Dirichlet boundary condition at the horizon and the Neumann boundary condition 
at infinity. 



7(±)| 



\r=r+ 



0, 



5.o(±)| 



0. 



(36) 
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that guarantee that the boundary conditions discussed above for $(r) are satisfied. In these equations 
the choice of q^~^' apphes if we want to study the case where $'"•* = in (33), while choosing the 



function q^ ' ahows us to address the case $(+) = 0. For the numerical implementation, it is convenient 
to use the variable 

y = l-^, (37) 

instead of the radial coordinate r, since y is dimensionless and bounded, < y < 1. We will also use 
a scaling symmetry of the radial equation to normalize all our quantities in units of the cosmological 



length. Then equation (29) or (34) depends only on four dimensionless parameters, namely, k, r^/i, 
fln^ and A± (or, equivalently, fxi). The strategy is now to fix k and run the spectral numerical code 
for several values of r^/£ and A± and find the eigenvalues Q.'j^i'^ of (34). These will be presented in 
Section [42;i 

As a double check, we will also use a shooting method to reproduce some of our results. A further 
motivation to use this second numerical method is that it allows us to find the important tail on the 
lower r_|_ region of Figure [5] with better accuracy. 

In the shooting method, the numerical strategy is the following. Again we want to fix k, r^/i 



and A± and find ^h^ at the onset of the instability. Equation (29) has two critical points: at the 
horizon and at infinity. We first focus our attention at the horizon. Using a Taylor series expansion, 
we construct the solution in the near-horizon region up to the eighth order in the radial distance to the 
horizon. Fixing k, this solution depends on r+/£, A±,Qh( but also on an arbitrary amplitude Ah- We 
then integrate numerically the radial second order ODE, using a standard fourth order Runge-Kutta 
method, up to a large radial distance. We repeat the procedure, this time at the asymptotic critical 
point where we start by obtaining the asymptotic solution up to eighth order. Again this solution is a 
function of r+/i, A±, 0,h£ and of an arbitrary amplitude A^o that we set equal to 1. We integrate this 
solution down to very small values of the radial distance. In the overlapping region of the two solutions 
we then do their matching. The requirement that both the scalar field and its radial derivative must 
be continuous fixes the values of the two unknowns, namely the desired O// and the amplitude Ah. 
The whole process is now repeated for a grid of values {r+/£ , A±} at fixed k to generate the results 
in Section [4X2] 

A subtlety intrinsic to the system at hand and not so standard in shooting applications is that the 
onset of our instability occurs close to the extreme solution. The Taylor expansion of the near-horizon 
solution is in terms of powers of y/ [Q'j^ — (fi^*)^] , where y is the compact radial coordinate defined 



in (37). Therefore, since the instability we search for sets in very close to extremality, typically the 
near-horizon expansion breaks down. To avoid this we work with a new compact radial coordinate 
defined as 



r2 



y = [ni-imrr{i--^)- (38) 



4.2.2 Results 



We expect an instability if we satisfy ( 32 ) . The simplest equal spin MP system that captures all the 
features of the scalar condensation instability is the d = 5 case with k = 0. Therefore we will present 
most of our results for this case, and in the end we will discuss the d >7 and/or k ^ cases. 

In d = 5, and k = one has ^^|sf^^ = —4, h'^Inhbf^'^ = —(3 + ^^/r^) and ^^|nrat^^ = —3, 
with ^j?\nh BF < fJ-'^\unit- Therefore both normalizable modes with fall-off r~ =■= should be unstable 



for masses obeying the condition (32) (for d>7 this is not true for large enough r+). This is precisely 
what we find numerically, following the strategy outlined in the previous section, as shown in Figure 
[l] In this figure, we plot the dimensionless angular velocity CIh^ (measured wrt the extreme angular 
velocity Q^^£) at which the HHT black hole has the zero-mode for the scalar condensation instability, 
as a function of the dimensionless horizon radius r+/£ and of the conformal dimension (s) A = A± 
defined in ([2]). The red curves describe the modes with faster fall-off r^ +, the blue curves describe the 
modes with fall-off r~ , and the single black curve is precisely at the BE bound where A_|_ = A_ = 2. 
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Figure 5: Detail of the black curve in Figure [T] where we fix the scalar mass to be at the BF bound 
i.e. when A+ = A_. The plot shows the dimensionless angular velocity wrt the extreme value, 
6M = (J^^* — ^h)^, as a function of the horizon radius, r^/i. 

The expectation is that for fixed A, all black holes that are in the parameter space below this surface 
are unstable and decay into hairy black holes with a scalar condensate surrounding the horizon. That 
is, the HHT black hole is stable for sufficiently high temperatures, but, for given r_|_ and A, there is a 
critical temperature where it becomes marginally unstable against scalar condensation, and for lower 
values of the temperature the black hole is unstable. Note that the surface extends to arbitrarily large 

r+- 

For clarity, we also plot the curve (il|^* — ^h)^ as a function of r^/£ for a fixed value of the scalar 
mass in Figure [5J for concreteness we picked the black curve in Figure [T] where A = 2. Note that the 
"tail" on this curve extends all the way down to an extremal solution with r^^/i = 1 (something that 
is not apparent from Figure [T|. But, for A = 2, this is precisely the threshold of instability predicted 
by the near-horizon BF bound. In more detail, for A = 2 the near-horizon BF bound is violated if, 
and only if, r^/i > 1. Hence we have numerical evidence that this bound is sharp: an instability 
appears as soon as it is violatedjj 

The instability against the scalar condensation is not the only instability present. Indeed, it 
competes with the well-known superradiant instability that affiicts rotating AdS black holes whenever 
Of/£ > 1 [I3FI Since there is competition between the two instabilities it is relevant to ask whether 
the scalar condensation instability occurs when il-n^ < 1- We find a negative answer. Recall that 
extreme black holes always have ^h^ > 1- Figure [T] shows that condensation occurs only when Qh is 
very close to ri|^'*, and we have checked that this implies that il//^ > 1. So the scalar condensation 
instability studied here always co-exists with the superradiant instability. 

We can now discuss what happens when we turn on the CP^ dependence defined in terms of the 
non-negative integer k. We find zero-modes of the instability for k > 0, at least up to k = 3. As the 
simplest example, in the left plot of Figure [6| we show the properties of the instability for d = 5 and 
At = 1. The stationary axisymmetric mode is indeed still present but the instability sets in for lower 
temperatures. Indeed the maximum difference between the extreme value Qfff of the HHT angular 
velocity and the angular velocity VIh where the instability switches on is smaller in the left plot of 
Figure [6] than in the k = case shown in Figure [TJ For k 7^ it is still true that the instability 

^ Just as for the hyperbolic black hole, the extreme solution does not admit a regular time-independent solution 
associated to the threshold of instability. If one takes the limit of the time-independent solution associated to non- 
extreme solutions then the value of the scalar field at the horizon diverges at extremality. 

^The superradiant instability is associated to perturbations that break the symmetry generated by d/dip. This 
instability is not confined to the scalar field, e.g., it is also present for gravitational perturbations. 
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Figure 6: a) Similar to Figure IT] but this time for modes with angular eigenvalue k = 1 (in the HHT 
background), b) Similar to Figure [I] but this time for the d = 7 MP-AdS background (and k = 
modes) . 



persists all the way up to the NH BF bound. To check this we take again as a reference point the 
curve A_|_ = A_. In this case, according to the condition (32), the zero-mode of the instability should 
kick in at extremality at the critical value of r+ given by (31). This value increases as k 7^ grows. 
For Av = 1, it is r^/i = 3. Our numerical results indeed confirm that the zero-mode of the instability 
is present for values of r_|_ at and above this value. This can be inferred by the range of r_|_/£ in 
the left panel of Figure [6], although a precise check requires analyzing the tail that develops close to 
extremality like we did for k = in Figure [5] 

Finally, let us discuss what happens as the dimension parameter d = 2N + 3 grows. We have 
explicitly checked the cases d = 7, 9 and the instability is still present. The numerical results for d = 7 
and K = are presented in the right panel of Figure [6] As d grows, the instability sets in for black 
holes that are closer to extremality, i.e. that have lower values of (fi^* — 0,h)/^^^- For k = 0, the 
instability persists all the way down to r^/i = 1 when A+ = A_, in agreement with the conditions 



(31) and (32) for the existence of the instability. 



4.3 Numerical results: nonlinear 

In the previous subsection, we found the properties of the onset of the scalar condensation instability 
in the codimension-1 Myers-Perry-AdS black hole. In this subsection, we include the backreaction in 
the gravitational field and construct the hairy codimension-1 rotating black holes that are associated 
with this instability. 



4.3.1 Hairy black hole ansatz. Equations of motion 

We want to look for hairy black hole solutions that reduce to the codimension-1 Myers- Perry black 



hole (24) when the scalar condensate vanishes, and that preserve the isometrics of the latter solution. 



Therefore we take the following ansatz for the gravitational and scalar fields, 

,2 



ds 



h{f 
<|) = <I)(r) , 



■^^'^idt^ + ^/rj^r, + fh{r)[di; + Aadx'' - n{r)dtf + fgatdx^dx^ , 



4f/(f 



(39) 



where, as before, Qab is the Fubini-Study metric on CP , and A = Aadx"' is the Kahler potential of 
CP . When the condensate vanishes, <I>(f) = 0, and g{r) = 1, (39) describes the MP solution (24) 
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after performing the radial coordinate transformation f ^^ r"^. This ansatz solves the equations of 
motion derived from the Einstein-scalar action ( 15 ) when the following equations are satisfied (taking 
A'^ = 1, i.e., d = 5 henceforth) 



{System of three coupled ODEs of 2'^'^ order for (/, h, <!>)} , 



(40a) 



2f3^2 ^f (^f^^y ^ ^2;j/2 ^ j^2 (2f2$/2 - 3)) - /if/' (3/l + fh')] 
/l2f3 [2(/i _ 4)^2 + f (^2^2$2 _ 12)] + 2£2C72 



(40b) 



n' 



-a 






(40c) 



where the system of three ODEs described in (40a) is cumbersome and thus we leave its explicit 
expressions for equation (91) of Appendix [BJ 

The construction of the hairy black hole amounts to determining f{r), h{f) and <^(f) that solve 
the system (40), i.e. (91), and (40c). Once these are known, g{r) is straightforwardly obtained from 



(40b). The constant C^ in (40c) is linearly proportional to the angular momentum of the solution, as 



we shall confirm later. 

The boundary conditions are the following. The function /(f) must vanish at the horizon, and we 
take this condition as our definition for the location of the black hole horizon. We require the hairy 
black hole to be asymptotically AdS and thus /(f) must approach 1 + f/l"^ asymptotically. For the 
same reason, h{f) must go to the unit value at infinity, while il(f) vanishes there. At the horizon 
we require 5, f^, h,^ to be regular. The scalar field asymptotic boundary condition is determined by 
the requirement of normalisability at infinity. Hairy black holes exist for scalar masses above the BF 
bound and below the NH BF bound. For concreteness, we work with the scalar mass /i = ^i\bf for 
which the numerics considerably simplifies. Summarizing, the boundary conditions for the non-linear 
problem are: 






0. 



0(1) 



/1-.00 ^ ^ + 1 + ^(^~"') ' ^l.=.i = ^(1) ' ^l.^oo - 1 + or') ; (41) 



a,""^ 



Ir— 5>oo 2f2 



7-~^\. 



$ 



0(1) 



(J) ~ 



+ 0{f~^). 



4.3.2 Results 

In this subsection, we present the scalar condensate, the temperature, the angular velocity and the 
entropy of the rotating hairy black hole as a function of either its energy or angular momentum. 
To construct the exact hairy black hole, we cannot use spectral relaxation methods, in contrast to 



what we did in Section 3.3 Relaxation methods crucially hinge on the positivity of the discretisation 
matrix representing the differential system at hand. When rotation is included, this is no longer 
possible, and the method at its most basic form does not converge. 



However, because (40a) is a one-dimensional system of non-linear differential equations, we can 



resort to a shooting method. Here we regard as the boundaries of our integration domain the horizon 
and the AdS spatial infinity. This method attempts to solve a boundary value problem, by reducing 
it to two initial value problems starting at each boundary. To determine the initial data at each 
boundary, we Taylor expand the equations in the neighbourhood of both the horizon and the AdS 
spatial infinity, determining the asymptotic solutions in a series expansion at each of the boundaries. 
We then integrate from both boundaries to an interior point in the integration domain, using a standard 
fourth order Runge-Kutta method, and demand that the two solutions match. 

A couple of comments, regarding the specific solution we are looking at, are now in order. Unlike 
previous non- linear hairy black hole solutions, such as the ones found in [2j and the ones studied in 



section 3.3, the rotating non- linear hairy black holes that we are going to study here only exist for very 
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small temperatures (Tf/£ ~ 10~^) - they are ultra-cold black holes. Furthermore, small variations of 
the temperature often correspond to large gradients of the entropy, as we shall see later in the phase 
diagram. Physically, the almost "infinite" throat characteristic of near-extreme black holes increases 
the proper length between the horizon and the cut-off scale where we choose to truncate the AdS 
space. This means that, in order to stabilise any numerical approach, we need a large resolution 
in both our discretisation scheme and moduli space of solutions. Finally, because the number of 
integration constants of the above system of equations is seven (recall that g is obtained algebraically 
once /, h and $ are known, and that Q obeys to a first order differential equation), the total number 
of parameters to be determined in the shooting process is seven. This makes the possible phase space 
of solutions very large, and thus very difficult to explore. Of course, these problems would be solved 
by using relaxation methods, but, as we have explained above, these are not available when rotation 
is included. 

To compute the energy of the hairy black hole in AdS we use the Astekhar-Das formalism |20j . 
The temperature, entropy, angular momentum and energy of these solutions are given by 



2.^h{rl)g{rl) ^ 4 \ & J 

(42) 
where 71 £^ is the 0{f ^) coefficient in the large r expansion of h{'r) and 72^^ is the 0{r ^) coefficient 
in the large f expansion of /(r). 

Ideally, we would like to present three-dimensional plots for the variation of the physical pa- 
rameters, such as the entropy and temperature, as functions of the energy and angular momentum. 
However, due to the difficulties alluded above, this does not seem feasible with the numerical methods 
we implemented. As such, we decided to determine the solution phase space at either lines of constant 
energy or lines of constant angular momentum. At the end of this section, we will explore the phase 
diagram along lines of constant horizon size r_|_ . 

For the HHT black hole, the energy along lines of constant angular momentum is bounded below 
by extreme solutions. Alternatively, the angular momentum along lines of constant energy is bounded 
above by extremality. We expect to see this zero temperature bound when analysing the phase 
diagram. However, we do not expect to find regular extreme hairy black hole solutions in the zero 
temperature limit. The argument is the same as we used for the hairy hyperbolic black hole (following 
Ref. |21|). An extreme black hole has a near-horizon limit. The scalar field must be constant in the 
near-horizon geometry, but then its equation of motion (for non-zero mass) implies that it vanishes 
there. Hence the scalar must vanish at the horizon of the full black hole solution. But we expect the 
condensate to grow, not decrease, as we lower the temperature. Hence the zero temperature solution 
cannot be a black hole. 

In Figs, hi we plot several physical quantities as a function of the dimensionless energy E/£'^ of 
the 5d rotating hairy black hole, for C^ = 3320 (i.e. J = 8307r^'^). In each of these plots there are two 
coexisting curves for some range of the energy. The range of energy considered is very small. This is 
because the temperature at which the bifurcation to the hairy solution occurs is so small that (from 
the first law) a tiny change in energy corresponds to a large change in entropy. The black curves are 
the numerical solutions corresponding to the rotating hairy black hole and the red curves correspond 
to the HHT black hole with the same angular momentum. There is excellent agreement between 
the linear results above and the non-linear solution when the condensate is sufficiently small, and 
in particular the point where the condensate vanishes coincides with the linear result with a 0.05% 
precision. To control the numerics, we have explicitly checked that the first law of thermodynamics is 
readily satisfied to a precision of 0.09%. 

The results are very similar to the hairy hyperbolic black hole: the condensate (or vev of dual CFT 
operator) becomes larger as the energy decreases, and both the temperature and entropy decrease with 
decreasing energy (in accordance with the first law of thermodynamics). When the hairy black hole 
has the same mass and angular momentum as a HHT black hole, it is the former that has the larger 
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(a) Condensate as a function of the energy. (b) Temperature as a function of the energy. 
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(c) Angular velocity as a function of the energy. (d) Entropy as a function of the energy. 

Figure 7: Different physical quantities for 5d rotating black holes as a function of the energy for 
J = SSOvr.^^. The black curve corresponds to our hairy black hole, the red curve to the HHT solution. 
The inset figures provide an expanded view of the region near the bifurcation point where the two 
solutions merge. 



entropy, which suggests that it should be more stable. However, from Fig. 7(c)[ we see that the hairy 
black hole angular velocity is always above £^^, indicating that it probably will suffer a superradiant 
instability [13j . Our code is numerically unstable for temperatures below 0.0004 .f~^, and that is the 
reason why we see an artificial lower bound on the energy before we reach the zero temperature 
solution. In analogy with the hairy hyperbolic black hole of the previous section, we strongly believe 
that the zero temperature and zero entropy solutions will coincide. 

In Figs. |8]we plot several physical quantities as a function of the dimensionless angular momentum 
J/£^ of the rotating hairy black hole, for E = 2717.85044.^^. (As mentioned above, the temperature 
is so small that a tiny change in the energy results in a large change in the solution, which is why 
we have to specify E so precisely.) Again, the smallness of the temperature implies that a small 
change in angular momentum results in a large change in the entropy, so the range of variation of 
J is small. Along the lines of constant energy, as the angular momentum increases, the condensate 
and angular velocity increase but the temperature and entropy decrease. The code is numerically 
unstable whenever we reach a temperature of the order of 0.0004.^"^, and this is why we never reach 
zero temperature. Again we observe that the angular velocity never decreases below i~^, and as such, 
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(c) Angular velocity as a function of the angular (d) Entropy as a function of the angular mo- 
momentum, mentum. 

Figure 8: Different pliysical quantities or 5d rotating black holes as a function of the angular momen- 
tum for E = 2717.85044 £^. The black curve corresponds to our hairy black hole, the red curve to the 
HHT solution. The inset figures provide an expanded view of the region near the bifurcation point 
where the two solutions merge. 

we expect these black holes to be superradiant-unstable. 

We have also probed other directions of the phase space, namely along lines of constant horizon 
size r+/L This direction explores a larger region of the phase diagram, i.e. both the energy and 
angular momentum vary more than in the previous directions, before reaching either the HHT black 
hole or the numerical cut-off temperature 0.0004 ^~^. In this way, we hoped to find black holes with 
angular velocity below £~^. However, that was never the case. In all directions that we probed, the 
angular velocity was always bigger than £~^. 

5 Rotating black holes: d = 4 Kerr-AdS 

At this point, it will not surprise the reader to learn that a massive scalar field can also condensate 
in the vicinity of a 4d Kerr-AdS black hole that is sufficiently close to extremality. In this section, we 
will determine the threshold mode associated to the condensation in this background. 
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5.1 Introduction 

The Kerr-AdS geometry is described by the hne element 

ds^ = -^ (dt -- sin' 9 dA\^dr^ + ^de' + ^ sin' efadt-'^^^dA , (43) 

S^ V ^ y Ar Aq S^ \ r. J 

where 

Ar = {r' + a')(l + '^)-2Mr, ^ = 1-^, A, = 1 - ^ cos^ , J:' = r' + a' cos' 9 . 

(44) 
This solution satisfies R^i, = —3i 'g^u, and asymptotically approaches AdS space with radius of 
curvature i. The ADM mass and angular momentum of the black hole are M/r? and J = Ma/r?, 
respectively [M|- The event horizon is located at r = r_|_ (the largest real root of A^). The angular 
velocity measured with respect to a non-rotating frame at infinity is 

The rotation parameter is bounded by 

a<L (46) 

Solutions saturating this bound do not describe black holes. In the limit a — t- ^ at fixed r_|_, the mass 
and angular momentum of the black hole diverge. The circumference of the black hole as measured 
at the equator becomes infinitely large in this limit. 
The temperature is given by 

r J- / rl \ 1 1 / r' 



^ 2^V^ l-^) rl+a? 47rrM^ "'^ ' ' ^^'''^ 



The extreme solution is the configuration with a = Oext; where 




«ext = r+^/:^^^ => nt = ^ :^, -, and ^<^. (48) 



Note that only "small" black holes with r+/£ < 3 ^'' can reach zero temperature in virtue of (46). 



The strategy is now similar to the one carried out in the HHT case. We study the Klein-Gordon 



equation for a massive scalar field in (43). We seek a stationary axisymmetric solution, which we 
expect to arise at the threshold of the instability. The explicit details of this study are given in 
Appendix [A} In the next subsection we present the results. 

5.2 Numerical results: linear 

The properties of the zero-mode of the scalar condensation instability are summarized in Figures [9] 



and 10 In the left panel of Figure [9j we plot again the dimensionless angular velocity ^h^ where a 
time- independent threshold mode appears (measured wrt Qff^^t), as a function of the dimensionless 
horizon radius r^/l and of the conformal dimension A defined in ([2|. The red curves (which combine 
to form a surface) describe the modes with faster fall-off r~ +, the blue curves (surface) describe the 
modes with fall-off r~ , and the black curves are precisely at the BF bound, which corresponds to 
A = 3/2. Note that, for low temperatures, (5J7 = {ytn — ^|/*) and r+ uniquely specify the solution 
(see Appendix [A]) . 

In this figure, there are two distinct surfaces present. These correspond to threshold modes with 
different angular dependence, analagous to the modes with different k that we discussed above. The 
scalar wave equation is separable in Kerr-AdS and one can label the solution by an integer / > equal 
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Figure 9: a) Threshold mode of the scalar condensation in the Kerr-AdS black hole. This plot shows 
the location of the threshold mode as a function of 5il.l = {iln — ^|f*)^; the parameter r^/i, and the 
conformal dimension A. b) Detail of the black curves in Figure a) with A = 3/2. The upper (lower) 
curve describes the I = {I = 1) harmonic. The dashed line corresponds to a = ^, i.e, infinitely large 
black holes. An instability with given / is present for points lying below the surface/curve with that 
value of /. 
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Figure 10: The harmonic structure of the zero-modes / = (left) and / = 1 (left) discussed in Figure 

El 

to the number of zeros of the angular part of the solution. The upper surface in the figure is for / = 0, 
the lower surface for / = 1. We find no numerical evidence for zero- modes with higher harmonics, 
/ > 2. Scalar condensation with given / occurs for points below the appropriate surface. 

In the right panel of Figure^ we plot in isolation the two black curves (i.e., with A = 3/2) of the 
left figure to see the details more clearly. This figure reveals that the instability occurs only above a 
minimum value of r+ corresponding to an extreme black hole. The curves end on a curve with a = £, 
corresponding to infinitely large black holes. 

In Figure 10 , we plot the amplitude of the I = scalar mode (left panel) and of the / = 1 scalar 
mode (right panel) as a function of angular and radial coordinates, x = cos 6 and y = 1 — r+/r. We 
clearly see that the / = harmonic has indeed no zero in the interval — 1 < x < 1 while the I = 1 
harmonic has precisely a single zero in the same interval. 

In the examples discussed above, we have seen that, at extremality, the near- horizon geometry has 
an AdS2 factor whose associated BF bound provides a sharp criterion for whether or not an instability 
is present. For Kerr-AdS, this is less clear because the near-horizon geometry is inhomogeneous. The 
radius of the AdS2 varies according to position on the S^ horizon. Hence there is no unique BF bound 
associated to the near-horizon geometry. 



22 



Finally we note that ^h^ > 1 whenever scalar condensation occurs so this instability coexists with 
the superradiant instability discussed in Refs. p9]-[33]. 

6 Stability results 

6.1 Introduction 

We have discussed several examples in which a free scalar field with negative /U^ is unstable in a 
near-extreme AdS black hole background. The aim of this section is to obtain sufficient conditions for 
such a scalar field to be stable. More precisely, given a black hole background we wish to determine 
a minimum value of fj? < which guarantees that the scalar field is stable. Our arguments will 
be based on those of Holzegel [TT], who proved stability of a free scalar field in the d = 4 spherical 
Schwarzschild-AdS background provided that the BF bound is respected, and (if below the unitarity 
bound) using boundary conditions defined by the fall-off r~ + . The main idea is to construct an 
energy functional for the scalar field that is non-negative and non-increasing. Hence if initially small 
then it must remain small. This can be used to prove decay of the field. We shall generalize only 
the first step of this argument (construction of a suitable energy functional) to various different black 
holes, and also to charged scalar fields. We expect that the second part of the argument (proof of 
decay) also can be generalized. 

6.2 Static black holes, uncharged scalar 

We shall follow the same approach as Holzegel but allowing for planar or hyperbolic spatial sections 
and general d. Consider a static, asymptotically AdS solution in d-dimensions, with metric of the 
form (It]) for some f{r). This encompasses both the Schwarzschild-AdS and Reissner-Nordstrom-AdS 
solutions. 

We define a new time coordinate t* by [11] 



dt* = dt+{j-Y^^^,]dr. (49) 



The line element ([7]) then reads 



ds' = -f dtl + 2 ^X'J^^i^J dUdr + ^^^^ij^i^^^l dr' + r' dT?, . (50) 

These coordinates are regular on, and outside the future horizon. They have the nice property that 
surfaces of constant t* are spacelike and intersect the horizon, with t* increasing along generators of 
the horizon. The Killing vector field timelike outside the horizon is 

Consider the energy of an uncharged scalar field on a constant time slice T^^ , with normal n'^ and an 
inner boundary at the horizon. The energy is given b}jj 

E{U)= [ T^.^n'^dS (52) 

where dS is the volume element on S^^ , and the energy- momentum tensor is 

T^, = {d^^){d,<^) -lg^,{{d<^f + f,''l>^) . (53) 



® In this section, we shall consider only boundary conditions corresponding to the decay r '^+ . If below the unitarity 
bound then one could also consider decay as r~ ^ . This requires modifying the definition of E. 
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On the background ( 50 ) , the energy is exphcitly given by 



dT,i 



r+ 



2 + 2r /e / ^2 ^ ^^,2 ^ J_(VcI,)2 + ^2^2 \ ^d-2 



(l+r2/£2)S 



dr . 



(54) 



where a dot denotes a derivative with respect to t,,,, a dash denotes an r-derivative and V is the 
connection on S^. 

Consider two constant time shoes at times 1^,2 > t*i- The bomidary conditions on the scalar field 
imply that the energy fiux at infinity vanishes, so conservation of the energy current implies that any 
change in the energy must result from a fiux across the horizon: 



E{t^,2) — E{t.ti) 



t*2 

t,i 



T^ueer''-^ 



dtifdYjji 



t.2 



$2 r'^-^ 



r=r^ 



dtjfdTik ) 



(55) 



where we used the fact that ^ is both tangent and normal to the horizon. The RHS is non-positive 
hence the scalar field energy outside the black hole is non-increasing. 

If /i^ > then E{t^:) is manifestly non-negative. If it is initially small then, since it is non- 
increasing, it must remain small ^^ Since E{t^:) is a sum of squares of $ and its derivatives, it follows 
that $ must remain small, hence there cannot be any scalar condensation instability. 

The interesting case, however, is //^ < with fj,"^ > /i2|^^. In this case, -E(t*) is not manifestly 
positive (the dominant energy condition is violated). However, we can exploit an argument of Ref. 
|llj to demonstrate positivity. We integrate by parts the last term: 



(j)2^d 2^^ 



r+ 



1 

d-1 



^2/^rf-l 



„d-l 



r+ 



2 
d-1 



$$'(r' 



W^d-l 



„d-l 



)dr . 



The surface term vanishes if we assume $ decays sufficiently fast at infinity, namely as r "• 
integral Hardy inequalitjrH implies that |llj 






„d-l 



dr . 



Hence, if ^u^ < we have 



where 



F{r) = fr''-^ + 



1 r°° 

E{U) >- F{r)<^'^drdJ: 

4fi 



d-2 



,2 



(d-1) 



^-id-2) fd^l 



„d-l 



(56) 
The 

(57) 

(58) 
(59) 



So if we can prove that F{r) is non-negative then the energy is non-negative. We can then argue as 
before: E{t^:) is non-negative and non-increasing and so must remain small if initially small. We can 



then deduce from (58) that $' must remain small and then (57) implies that ^ must remain small so 
there cannot be a scalar condensation instability]^^! 



We now examine the form of F{r) for Schwarzschild-AdS. Using r > r^ and /x < gives 

,2 



F{r)/r > 1 



/^ 



u2 



,d-l 



rt^ \+k{r 



d-3 _ d-3 



(60) 



Hence if jj? > fJ-'^lnp then F{r) is manifestly non-negative for k = 0,1. (This generalizes the argument 
of Ref. [llj to general d and to A; = 0.) Hence a scalar field in the spherical or planar Schwarzschild-AdS 
background is stable if it obeys the BF bound. 



^"This is a standard argument which can be apphed whenever the dominant energy condition is satisfied. 
^^This can be derived using the Schwarz inequality on the RHS of (561. 
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Ref. [TT] gives more rigorous arguments to show boundedness of $ 
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For the k = —1 case, a positive integrand in (58) requires 



M 



2 



^1 > max < :r-i — }■ (61) 

Sketching the function in curly brackets on the RHS of this relation for a non-extreme solution, it 
starts at — cxd at r = r-|_ then increases to a negative maximum value greater than —1, then decreases 
monotonically to — 1 as r — t- oo. Hence F{r) > implies that ^^ must lie strictly above the BF bound. 
(For large black holes r_|_ > 2i (say), the value at the maximum is only slightly greater than —1 so 
the new bound is very close to the BF bound.) However, as r^ — )• r'^^, the location of the maximum 
tends towards r = r.^- and the above inequality reduces to 

2\ «— 1 2\iSchw) .„„^ 

^ -"l£^=^ \nhbf- (62) 

In other words, F{r) is non-negative if, and only if, the scalar field respects the BF bound associated 
to the near-horizon AdS2- Hence if this bound is respected then there can be no scalar condensation. 
This argument demonstrates that the extreme black hole is stable if the near-horizon BF bound 
is respected. But is this bound sharp? Maybe, although it is not manifest, the energy is still positive 
for even lower //^. This is easy to exclude: we shall exhibit a trial function $ for which the energy is 
negative if the near-horizon BF bound is violated. For a general black hole with metric of the form 
([7| extremality implies that 

/(r) = (r - r+f/R^ + 0((r - r+f) . (63) 

where R is the AdS2 radius. Consider a trial function ^eif) defined by 

$e(r-) = Ce(r-r+ + e^)-^/V^, (64) 

where C^ is a normalization constant fixed by demanding that Vs f dr <1>^ fd,-2 _ j^ ^y^ ^g ^j^g volume 
of the unit radius compactified hyperboloid) . This trial function is regular at the horizon and satisfies 
the decay conditions at infinity. For small enough e, the associated energy is dominated by the 
near-horizon contribution, 

E, = l{f^'- f^'lNHBF) + O (y^^^) • (65) 



Therefore, for any extreme black hole of the form (I7|), if /i < fi \nh bf, then there are regular initial 
data for which the energy functional is negative. For the hyperbolic Schwarzschild-AdS black hole, 
our numerical results reveal that an instability is present in this case. 

6.3 Rotating black holes 



We focus now on the rotating black hole solutions (24). Let us introduce a coordinate system analogous 



to ( 49 ) which is well behaved on the horizon. Consider new coordinates t* and ip^ , such that 



dU = dt+{j-Y^^jj^j^dr, dip, = dip + i^-j- -^^^2112 ]^hndr. (66) 



The line element (24) then reads 



,2 / , 2 2 1+ r^/f - / , , 2 + 2rV£2 _ / 
ds^ = - f dtt + —F- —. r, ,„^ dUdr + — '„r,J dr^ 



h * ^/h l + r^/P * (l+r2/£2)2 

+ r'^hldip^ + Aadx"" - Qduf + r'^gabdx^dx^ . (67) 
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The Killing vector C = 9i, + Vtnd^^ is everywhere timelike outside the horizon for |r2//|^ < 1. Let 
us define the energy on a constant time slice St, with normal n^ as E = J^ Tfj^jj^^^n^ , i.e. we follow 
Ref. [29] and define the "energy" with respect to the Killing vector ^, as opposed to dt^ 
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On the 



background (67), the energy is then explicitly given by 

"2 + 2rV£2_y 



E 



-TT 



fiS 



N 



r+ 



(l+r2/£2)S 



h{c<^y + 



1 



(0^-0)2 (9^.$) 



+ f{dr^ + 



2 + 2rV^2 _ f f^ 
(l + r2/^2)2 j< 



2^2 



inH-n)d^,'^] +^(p'^y + ^^'^ 



r'^-'^dr . 



(68) 



where Va = Va — Aad^^, while V is the connection on CP , whose integration measure we denote 
as Etv- As before, it is straightforward to demonstrate that £^ is a non-increasing function of t=K (this 
uses the fact that ^ is tangent to the generators of the horizon) . 



If // > 0, every term in the integrand of Eq. (68) is non- negative when 1^1/^ |^ < 1. Hence there is 



no scalar field instability [IH]- When jf^iifl^ > 1, the coefficient in front of (9^,<1>)2 becomes negative 
in a region near infinity, where ^ is spacelike. This is the signal for the superradiance instability |13j . 
which affects perturbations breaking the d^^ = d^ rotational symmetry. 

In the present paper we are interested in axisymmetric perturbations, i.e., d^^^ = 0. In this case. 



the relations (58), (59) and (61) still hold (with dS^ substituted by 2TrdT;N)- Hence we have a lower 
bound (61) on ^"^ which guarantees non- negativity of E and hence stability against axisymmetric 
scalar perturbations. In the extreme limit, the maximum on the RHS of (61) is located at r = r+ for 



large black holes, and at spatial infinity for small black holes. The energy is non- negative for 



/^ 



> 



N + l 



N + 2 + N- 



2, (MP) 

/^ \nhbf 



if 



r+ > 



(69) 



The trial function argument (63)-(65) shows that this bound is sharp in the sense that there exist 



negative energy initial data if it is violated. Our numerical results confirm that an instability appears 
if this bound is violated. For small black holes, r_|_ < £, the inequality is simply ^2 > 



Z^^Ibf, 



I.e. 



these solutions are stable against axisymmetric scalar field perturbations as long as the asymptotic 
BF bound is respected. 

6.4 Charged scalar field 

Since the condensation of a charged scalar field in Reissner-Nordstrom-AdS has attracted so much 
attention recently, we shall consider the extension of the above arguments to this case. The extension 
is not completely straightforward because the charged scalar field current couples to the background 
gauge potential. The equation of motion for the scalar field is 



L>2$ 



/i'^, 



(70) 



with D^ = V^ — iqAf^, where q is the charge of the scalar field and A^ is the gauge potential. We 
consider here Reissner-Nordstrom black holes with line element ([7]), where 



f{r) =(k- 
and the gauge potential is given by 

A = 



(r r+)"'~'^ 



r.d—3 



+ 




rd—3 



dt 



Q 



jr 



d-3 



,.d-3£2 



dt, 



(71) 



(72) 



'^^ This means that our "energy" is actually £ — Q.hJ, where £ is the energy defined with respect to d/dt and J is 
the total scalar field angular momentum. 
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where 7 = yj2{d — 3)/(d — 2). We are using the gauge in which the potential vanishes at the horizon. 



The hnear instabihty for scalar condensation can be found by solving Eq. (70) leaving the gauge 
potential unperturbed [Ij, since this latter perturbation would have a higher order contribution. How- 
ever, energy considerations require a more careful treatment. As noticed by |35] . the charged current 

j^ = iq[{D^<^)^^-{D^^)<^^] (73) 

sources a gauge potential perturbation through 

y^F^, = -j^. (74) 

The consequence is that the energy-momentum tensor for the scalar field, 

T^u = {D^^){D,<^)^ + {D^^)\D,^) - g^u{\D^'' + ^x^\^^) , (75) 

where | • P = {•)\-)-, is not separately conserved, V^T^u = F^^j^^. Only the total energy-momentum 
tensor, which includes the contribution from the perturbed gauge potential, is conserved. Therefore, 
the energy current J^ = T^^u^'^, where ^ is the timelike Killing vector, is not conserved either, V^ J^ = 

This can be easily fixed by considering instead the energy current (in the Lorentz gauge V^A^ = 0) 



Jf, = T^ue-JMun, (76) 

which accounts for what would be the perturbed gauge potential contribution. Note that this is not 
gauge invariant. The conservation of this energy current is guaranteed by the condition A^"'^''' = 0. 



Using the coordinate system in (50), the energy E = Jj, J^n^ of the charged scalar field is given by 



E=id^,j (' ' + ^^'{^^,)/ \D^M' + f\DM' + ^,m? + ^'m' - /*^.,) r'^-'dr , (77) 



^■4 



whose only difference with respect to (54), apart from the fact that we have a complex scalar field 
and gauge covariant derivatives, is the presence of the last term. 

The energy fiux through the horizon is determined by J • ^ = Fiu^^^'^ which is easily seen to be 
positive at the horizon. Hence the energy is non-increasing in time ^^ 

The squares in the integrand can be rearranged to give 

> J d^^ £ (tttw^ '^*-^'' + ^'^^•^'' + ^'^^'' + ^'"''0 ^'"'^ ' ('') 

where ^g = iJ,'^ — max{/~-'^g^A^^}. 

If Mg — then the energy is manifestly non-negative. If //^ < 0, we can proceed as in the uncharged 
case. Notice first that \dr'^\^ > (9r|$|)^ Hence, if /x^ < 0, 

/•oo 

EiQ > / F,{r){drmfdrdEk, (79) 

Jr+ 



^* There is a close analogy with the rotating case here. The "standard" choice of gauge would be A^ = at infinity. 
This leads to the "usual" notion of energy £. This is analogous to the energy in the rotating case defined with respect 
to the usual generator of time translations d/dt. However, to prove stability we make the alternative gauge choice, for 
which Afi = at the horizon and we call the energy E. This corresponds to the energy in the rotating case defined with 
respect to the co-rotating Killing field ^. In the charged case we have E — £ — Q(f>, where Q is the total charge of the 
scalar field, which is analogous to the result E = £ — QhJ in the rotating case. 
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where 



F.ir) 



fr' 



.d-2 



+ 



^^g -(d-2) ( d-i _ d-i^ 2 
{d-lf V "^^ 



A positive integrand in ( 79 ) requires 



/^g 



tJ''^\BF 



> max 



fj^^2{d-2) 



(r 



d-l 



d-l\2 



<-) 



(80) 



(81) 



This gives a sufficient condition for stability against charged scalar condensation. 
Consider now extreme black holes. These solutions satisfy 






d-3 
d-l 



(7202 _ ^) . 



(82) 



The maximum of / q Af^ is located at r = r_|_, and thus /Xg is the near-horizon "effective mass", 



^g 



2j,2 



^i'-if-'q'Ai)r=r+=^' 



{d - 3)q^(l) 
2(d - 3)(/>2 - A: 



(83) 



The maximum on the RHS of (81) is also located at r = r^ if A; = or A; = —1. On the other hand, 



if A; = 1, that maximum is located at spatial infinity for small r^/£, otherwise it is located at r = rj^. 

d-3 



We get from (81) the inequality 



^2 > _d^ i^(^, _ ^^^2 _i.\^ „2|(fl^) 



4H 



(2((i - 3)(/>2 - A:) =fi^ 



\NHBF 



if A; = 0, -1, or if A; = 1 with ^ > ^ " , (84) 

t yd — 1 



which is analogous to the uncharged cases (with ^u — t- fig). This is shown to be a sharp bound by the 
trial function argument (63)-(65), i.e., there exist negative energy initial data if it is violated. 



Small spherical black holes are an exceptional case. In this case, (81) reduces to 



Mq >/^%F 



if A; = 1 with 



< 



d- 



(85) 



since in this case h'^\bf > A* 



2dRN) 
\NH BF- 



The difference to the uncharged cases is that this bound can be 



violated, without violating the BF bound for AdSd, because fig < /i^. Small extreme charged black 
holes can be unstable to the condensation of a charged scalar field even if the near-horizon "effective 
mass" Hg is above the near-horizon BF bound. Such an instability was found in Ref. [9]. The result 



(85) is unlikely to be sharp: it is sufficient but probably not necessary for stability. 

Let us make two final comments. First, notice that the charged scalar condensation on a charged 
black hole corresponds to the (charged) superradiance instability in AdS. The threshold of scalar 
condensation is a time-independent mode in the gauge (72), where A{r^) = and A(+oo) = 0. To 
relate to the known phenomenon of superradiance, we consider instead a gauge where the chemical 
potential vanishes at spatial infinity, A = A — cpdt. The scalar field transforms as <^ ~ ^-ruit _^ 
$ ~ e~*"*, where uj = u + q(j). The threshold of superradiance is w = q(p, which corresponds to a 
time-independent mode (w 

oo. In this case, only A; = 1 solutions 



0) in the gauge (72). 
The second comment is on the asymptotically flat limit i 



exist and fj? must be non- negative. A bound simpler than (81), but sufficient for the present purpose, 
is obtained for a scalar field satisfying fi^ > q'^j~'^, which corresponds to a BPS-like bound. In this 
case, the energy (78) is non-negative if / > 7^^t^- Since 



/ 



-7242 
7 ^t* 



fl 



1'^') 



^d-3 



i-3 



rd—3 



(86) 



the energy is non- negative if j'^cj)'^ < 1. The saturation of this bound corresponds to extremality, as 
seen in the £ — >• 00 limit of (82). Therefore, there is no instability of such a scalar field. 
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A Analysis of Kerr-AdS black holes 

The Kerr-AdS black hole has the peculiarity that, for small values of the dimensionless horizon radius 
r+/^, as the rotation parameter a or the ADM angular momentum J increases, the angular velocity 
Qh first increases and then decreases. Stated in other words, it is possible to have two different values 



of a with the same value of the angular velocity ^h- Recall that Qh is given by (45) and it is the 
angular velocity measured by an observer in a frame that does not rotate at infinity. We do not find a 
detailed discussion of this property in the literature. We discuss explicitly the Klein-Gordon equation 
for a massive scalar field in the Kerr-AdS background whose study generates the results presented in 
Section [H 



For r_|_/£ < 3 ^'"^^ as the rotation parameter a increases, the angular velocity Qh in (45) first 



increases, then reaches a maximum at a = r+ and then decreases monotonically until a = Oext- Along 
this path the temperature always decreases until it vanishes. For r^/£ = 3~^", one has 0|^* — )• 1/i 
when a — )• Oext = ^i and this is a singular zero temperature configuration. For 3~^'^ < r^/i < 1, 
Qh also has a maximum at a = r+ and then decreases to iln — ^ 1/^ as a — ?■ ^. Along this path the 
temperature decreases but never vanishes. Finally, for r^/i > 1, the angular velocity always increases 
monotonically and the temperature decreases without vanishing as a approaches i. This discussion is 
illustrated in Figures [Tl]|12[ We plot the dimensionless function i^n^ as a function of the dimensionless 
rotation parameter a/i for the five distinct cases that can describe a black hole as r^/i runs from 
small values to large ones. In all these plots, as a/i grows the angular momentum increases and the 



temperature (47) decreases. Because of the constraint a < I, only the black holes of Figure 11 (left) 



can reach extremality. 



We can discuss these properties in a complementary way. Using (45), we express the rotation 



parameter a in terms of the angular velocity Q^. There are two possible roots: 



r: 



2 I e2 



+ e"^ ± J {rl + £Y - "^^l^^^l 



For r^/£ < 1, one has a{Q}{) = O- while a < r_|_, and a{Q}{) = o+ for r^ < a < minjaextj^}- For 



r+/£ > 1, one always has a(0/f) = a-. In Figures 11 12 we also identify which of the roots a± 
describes the black hole in each situation. 

For reasons that will be soon clear, we will find it useful to work with a different quantity i^H 
defined in terms of 0,h by 



Qh = 73 V (^+ + ^^) ~ 4r^£^0|^ , with e = ±1 if a ^ r+ , respectively. 

This quantity decreases continuously and monotonically as a grows from zero to minjoext)-^}- It is 
positive if a < r-|_ and negative when r+ < a < minjacxt,^}- 



Figure 13 on the other hand summarizes some of the information of the previous plots: it gives a 
region plot of the parameter space where Kerr-AdS black holes exist. We express this information in 
terms of the variable iln that is relevant for our eigenvalue problem below. 

Consider the Klein-Gordon equation (jTl) for a massive scalar field in the Kerr-AdS background. We 
are interested only in the stationary and axisymmetric zero- mode of the scalar condensation instability, 
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Figure 11: Dimensionless quantities ^h^, ^h^, o,±/i as a function of the diniensionless rotation pa- 
rameter a/i for r^/£ < 3~^" (left) and r^/i = 3^^'^ (right). The transition between the curves 
a = a^ and a = a- occurs at a = r^. 




Figure 12:_Similar to Figure |ll| but for 3"^/^ < r+/i < 1 (left) and r+/£ > 1 (right). Note that for 
r+/e = 1, Q// = at a/i = 1. 



so we look for solutions $ = $(r, 0) ° Introducing the compact radial coordinate, y = 1 — r^/r, and 



the angular coordinate x = cosO, we find that the Klein-Gordon equation for ^{y,x) can be written 
as ^ 

L^{y,x) =nHiA^{y,x), 

where the second order differential operators L and A are given by (2/4. = r^/i) 



(89) 



L = 4(l-y)(l + y2) 



A = 4(1 -y) 



(1 - y)y [(1 - yfy + (4 - y{5 - 2y))yl] d^ 
-2 [(1 - yfy + yl{2 - y)[l - (1 - y)y]] dy - {l - x^) (1 - y) {l - ylx^) dl 
+2x(l -y)[l + yl (1 - 2x2)] a, + A^JLHW+ (l + vl) [l + x\l - yf] , 

y[2 - (3 - y)y] [(1 - yf + yl] dl + 2 [(1 - yf + y^] 9, 



(90) 



+ (1 - x^) (1 - y) (1 + yix^) dl - 2x(l - y) [l - y^ (l - 2^2)] 9, - 4^^^^^^ ^-^ _ ^2^-^ _ ^^2] _ 



^^It is well-known that, under the ansatz $(r, S) — R{r)S{6), the Klein-Gordon equation in Kerr-AdS can be further 
separated into a radial equation for R{r) and an angular equation for S{9). The latter has the AdS spin-0 spheroidal 
harmonics as solutions, S{9). Unfortunately, these harmonics cannot be constructed analytically, neither can the associ- 
ated angular eigenvalues of the separation constant. It is possible to construct these approximately in a series expansion 
for small values of the rotation and cosmological constant, but we are interested precisely in the opposite regime of large 
rotation. For this reason we do not find any advantage to introduce this separation ansatz. We search directly for the 
zero-modes of <l?(r, 6). 
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Figure 13: Region plot of the parameter space where Kerr-AdS black holes exist (area in between 
the several curves), expressed in terms of the eigenvalue variable ^u^ and the dimensionless horizon 
radius r+/£. Note that the / = and / = 1 zero-mode curves described in Figure [9] have as endpoints 
the curves labeled as 0+ = Ccxt and a+ = £ in the current figure. 



The ultimate reason why we introduced the quantity ^u instead of working with the black hole 
angular velocity ^u is now clear. Indeed, in terms of the former quantity, the Klein-Gordon equation 



(89) is explicitly an eigenvalue equation for ^u- The reason to follow this strategy is motivated by 



the availability of numerical techniques for solving eigenvalue equations of the form ( 89 ) , namely the 



spectral method already described briefiy in Section [4| We normalize all quantities in units of the 
AdS length. The Klein-Gordon equation (89) depends on the three dimensionless parameters r_|_/£. 



^H^ and A (or, equivalently, n€). Running a spectral numerical code for several values of r+/^ and 
A, we find the eigenvalues tint of ( |89[ ), and thus the associated physical angular velocity Vlui where 
the zero- mode of the instability appears. The results are presented in section [5} 

B Coupled system of ODEs for the codimension-1 rotating hairy 
black hole 

In this Appendix we write explicitly the system of three coupled 2nd order ODEs for {/, /i, <!>} that the 
codimension-1 rotating hairy black hole must satisfy and that are symbolically described in equation 
(40a[) of the main body of the text. These are: 



= {3/f2 [/i2f3 ^2{h - 4)f + f (/i2£2$2 _ 12)) + 2e^C^] } ^ 
efCl [Sfh + f {2fh' - f (3/i + rh') + /f (2/i$'2 + h"))] 
+ r^ hff (3/i + fh') (l2(/i - l)hf + f (6(/i - 2)f + r (/i^^^^^ _ ^2)) h') 

+ f[- ^hf"^ (6(/l - 2)^2 + f (;^2^2^2 _ ^2)) /j'2 j_ ^3 (^_g(^ _ 2)^2 _ ~ (^2^2^2 _ -^2)) /^/3 

- 2h^fh' [ - 6(/i - 4)£2 - 3f (^2^2$2 _ ;l2) + f2 (g(/j _ 2)£2 ^ ^ (^2^2^2 _ ^2)) $'2j 



+ 3/1^ (-4(/i - l)hf (2f2$'2 - 3) + f2 (2(/i - 4)£2 + r {^?f^'^ - 12)) h") ] ^ = , (91a) 

= 3/f [/i2f3 (2(/i - 4)^2 + ^ (^2^2^2 _ 12)) + 2fCl] f" 

+ f/'2 [-GfC^ + hr^ (l2(/i - l)hf + f (6(/i - 2)f + f {fx^f^^ - 12)) h')] 
+ /2f2 (_12(/l - 1)^2 ^ ~ (^2^2^2 _ -^2)) (_3/i2 ^ ~ (^f^> (^^ ^ ~^/^ ^ 2/iV$'2)) 

+ ff'UfCl + f^[2hf {3hf - r (^^£2^2 _ -^2)) h' + f2 (_6(_2 + h)f - f {fi'^f^'^ - 12)) h'^ 
+ h^ [l2(-8 + 5h)f + 3f (^2^2^2 _ 12) - 2f2 (6(-2 + h)f + f {fi^f^"^ - 12)) ^'^j ] | , (91b) 
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= 6/f [h^r^ (2(-4 + h)e + r {fiHH^ " 12)) + 2fCl] $" + 9fh^rH^f,^^ 



3£V^^ 



+ /ifV (3/1 + rh') [3^V^^ + 2 (6(-2 + h)f + f (/x^^^^a _ ^2)) $'] _ /j/if^/i' 

+ 2 (6(/l - 2)^2 + f (^2£2$2 _ ^2)) $' + f5;j/2 [3£2^2^ ^ ^ (6(/i - 2)^2 + f (^^^2^2 _ ^2)) $'] 

6^2c2 + /i2f 3 ^12(5 - 2/i)£2 + f (72 - 6£2^2^2^ 



+ 2$' 



+ f2$' (3^2^2^ ^ 2 (6(/l - 2)^2 + ^ (^2^2^2 _ ^2)) $') 
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